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(-■h ! Abstract 

With a method m which the Friedmann equation is written in a form such that evolution of 
the scale factor can be treated as that of a particle in a "potential" , we classify all possible cosmic 
^ij^ evolutions in the DGP braneworld scenario with the dark radiation term retained. By assuming 

that the energy component is pressureless matter, radiation or vacuum energy, respectively, we 
" find that in the matter or vacuum energy dominated case, the scale factor has a minimum value 

. . .e .... .e, .e „ . . 

circumstances, and there may exist an oscillating universe or a bouncing one. If the cosmic 
. scale factor is in the oscillating region initially, the universe may undergo an oscillation. After a 



(N 



X: 



number of oscillations, it may evolve to the bounce point through quantum tunneling and then 
expand. However, if the universe contracts initially from an infinite scale, it can turn around 
and then expand forever. In the vacuum energy dominated case, there exists a stable Einstein 
static state to avoid the big bang singularity. However, in certain circumstances in the matter 
or vacuum energy dominated case, a new kind of singularity may occur at ao as a result of the 
discontinuity of the scale factor. In the radiation dominated case, the universe may originate 
from the big bang singularity, but a bouncing universe which avoids this singularity is also 
possible. 
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I. INTRODUCTION 



The modified gravity has spurred an increasing deal of interest recently, because it can 
explain, without the introduction of an exotic dark energy, the present accelerating cos- 
mic expansion discovered firstly from the Type la Supernovae (Sne la) jl-6|. Among the 
modified gravity theories, the Dvali-Gabadaclze-Porrati (DGP) braneworld scenario 
generalized firstly to cosmology by Deffayet [8], is a very simple and popular one. The 
DGP theory starts with the idea that our observed four- dimensional Universe resides in 
a five-dimensional, infinite-volume Minkowski bulk and the whole energy-momentum is 
confined on a three dimensional spacial brane. In contrast to the Randall- Sundrum 
and Shtanov-Sahni lOj braneworld scenarios with high energy modifications to general 
relativity, the DGP brane produces a low energy modification (for a review of the phe- 



nomenology of the DGP model, see Ref. [11]). 

Since there are two different ways to embed the 4-dimensional brane universe into the 
5-dimensional spacetime, the DGP model has two separate branches denoted by e = ±1. 
The e = +1 branch is self- accelerating in the sense that the universe is rendered to 
accelerate at late times due to the lowly leaking of the gravity off our four-dimensional 
world into an extra dimension and to evolve eventually into a de Sitter phase jsj. However, 
the e = — 1 branch is very different since it does not self-accelerate. Thus, in order to 
explain the present cosmic acceleration in this branch, dark energy is required on the 



brane, like in the LDGP model 



12j and QDGP model (l3|. 



The inflation and preheating on the DGP brane have been discussed in Refs. 14j-ll9| 
and some new characteristics have been found. For example, the DGP inflation driven by 
a single scalar field with an exponential potential yields much better consistency with the 
current observation data |il5j|. Recently, we discussed the stability of the Einstein static 



universe in the DGP scenario 



20| and obtained that the universe can stay at this stable 



Einstein static state past-eternally, undergo a series of infinite, non-singular oscillations, 
and then evolve to inflation. Therefore, the big bang singularity can be avoided. Moreover, 



the cosmic background evolutions in the DGP model have been studied in 



[lO, 



21 



With a large value of the dark radiation term, it was found that the spatially flat DGP 
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braneworld gives the same dynamical possibilities of the cosmic evolution as a closed 
FRW universe and these possibilities include the oscillating, the bouncing, the Einstein 
static universes and the so-called loitering universe (see Fig. (4) of 22|). In the present 
paper, we plan to classify all possible cosmic evolutions in the DGP braneworld with a 
method in which the Friedmann equation is written in a form such that evolution of the 
scale factor can be treated as that of a particle in a "potential". The effect of the dark 



radiation are also considered in contrast to Ref. [20|]. Different from Ref. 2^, we keep, 
in our discussion, the spatial curvature term and do not impose the condition of a large 
value of the dark radiation. Let us note that this method has been used to classify the 
cosmic evolution in the Horava-Lifshitz gravity {23 1. 



II. THE FRIEDMANN EQUATION IN THE DGP BRANEWORLD 



We consider a homogeneous and isotropic universe described by the Friedmann 
Robertson- Walker (FRW) metric 



-dr + a\t) 



+ r^d^n 



(1) 



1 — kr"^ 

where a is the cosmic scale factor and k is the constant curvature of the three-space of the 



24| 



(2) 



FRW metric. In the DGP brane scenario, the Friedmann equation can be written as 

+ ^ = ^[P + Po(l + e^(p, a))] , 

where H is the Hubble parameter, p the total energy density and /i a parameter denoting 
the strength of the induced gravity on the brane. A is given by 



A 



Po V 



1/2 



(3) 



where 



^0 



with A defined as 



l-2r] 



p2A 
Po ' 



6ml 
PoP^ 



{0<T]<1] 



4 m^ 
PQ = m^ + 6^ 

P 



(4) 



(5) 



Here K5 is the 5- dimensional Newton constant, ^^^A the 5-dimensional cosmological con- 
stant in the bulk, A the brane tension, and Sq an integration constant related to the Weyl 
radiation (dark radiation) which is assumed to be positive in this paper. For simplicity, we 
restrict ourselves to the Randall- Sundrum critical case, i.e. A = 0, then Eq.(2) simphfies 
to 



In the very early era of the universe the total energy density should be very high. Thus, 
we will, in the following, only consider the ultra high energy limit, p ^ Po- In addition, 
we let 7] = 1. As a result, the Friedmann equation reduces to 



The dark radiation term is retained here in contrast to Ref [20i] where it is neglected. 
When a is small this term is very important. It is easy to see that the above equation 
describes a 4-dimensional gravity with minor corrections, which implies that p must have 
an energy scale as the Planck one in the DGP model. 

For the cosmic energy, we assume that it has a constant equation of state u and thus 
its density can be expressed as 

where g is a. constant. In the following, we take u = —1, 1/3 or 0, which corresponds 
to the vacuum energy, radiation, or pressureless matter dominated universe, respectively. 
Thus the Friedmann equation becomes 



^^ + ^-3^(^ + y2Po(^-.^§)). (9) 

Clearly, a > 

= (£^)i/(i-3.) is 

required when u 7^ 1/3, which means that the universe 



9 

begins to evolve at a > Oq rather than a = 0. So, the classical big bang singularity can be 



avoided. Let us note that this finite size initial universe can be created from "nothing'' 



through quantum tunneling 



25 



26| . For u = 1/3, g > p^£^o is needed and in this case 



the universe can originate from a = 0. 



Now we rewrite the Friedmann equation in the following form 

+ V[a) = 0, 



(10) 



where 



V[a) = k 



1 g e 




3a)— 1 
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(11) 



Thus V can be regarded as a "potential" and the scale factor a chang particle 
moving in it. This "potential" must satisfy the condition V{a) < 0. This gives the 
possible range of a when the universe evolves. Therefore, we can classify the types of the 
universe by the signs of k and e, and by the values of other parameters. 
All cosmic evolution types in the DGP braneworld are: 

(1) [Bounce]: If V{a) < for a e [qt, oo) and the equality holds at a = qt, a spacetime 
initially contracts from an infinite scale, and it eventually turns around at the finite scale 
ay, and then expands forever; 

(2) [Oscillation]: V{a) < for a G [cijnim (^max] and the equality occurs at a — otmin and 
o = dmax, a spacetime oscillates between two finite scale factors; 

(3) [FS oo]: V{a) < for a e [oq, oo). The universe starts at finite size(FS) oq and 
expands forever. 

(4) [BB =^ BC]: V{a) < for a e (0,0^] and the equality holds at a = ar- A 
spacetime starts from a big bang (BB) and expands. It turns around at a = and then 
contracts. Eventually, the universe contracts to a big crunch (BC). ot is the scale factor 
where the universe turns around from expansion to contraction. 

(5) [BB =^ oo or oo =^ BC]: V{a) < for any positive values of a, a spacetime starts 
from a big bang and expands forever, or the spacetime always contracts to a big crunch. 

(6) [FS =^ FS]: V{a) < for a G [ao,am] and the equality holds only at a = am- 
A spacetime starts from a finite scale ag and expands. It turns around at a = and 
begins to contract. When the universe contracts to the minimum scale a = ao, it should 
expand again. However, its evolution will be discontinuous at a = oq since the potential 
V{ao) 7^ 0. Thus, there exists a new singularity in this type. 
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III. THE EVOLUTION OF A MATTER-DOMINATED UNIVERSE IN THE 
DGP BRANEWORLD 



If the universe is dominated by pressureless matter [u = 0), the cosmic energy density 
can be expressed as p = Thus, the potential becomes 

V{a) = k - -i^^ - -i^V2po(^?™a-/i2£o) • (12) 

Clearly, oq = ^^-^ and H{aQ) ^ in general except for the case 

gl = 3kfi% . (13) 

This condition gives a boundary to obtain an oscillating universe. 

A static universe appears if there is a solution a = as > ao which satisfies V{as) = 
and V'{as) = 0. At as, both the cosmic expansion speed and acceleration equal to zero 
and thus the universe can stay at this point if it is stable. Differentiating V{a) with 
respect to a, we have 

V'{a) = - 1^9^^ (^4) 

Combining V{a) = and V'{a) = 0, we obtain, to get a static universe, a relation between 
gm and other parameters: 



9m = 9m=^^^ (^9A;/x(/i2 + 2£:opo) ± v^(3p2 - 28^p^)i^ , (15) 

which gives another two boundaries for obtaining the oscillating universe. Now, we have 
three boundary conditions (Eq. (fT3l flSil ) for an oscillation. Using Eq. (ITSll and Eq. ( fT2l) . 
one can find the static state solution 

as = % = , = , (16) 

^J9kfi{fi^ + 2£oPo) ± V3(3p2 - 2£oPo)^ 



which is a double root of the equation = under the condition V^' = 0. Then the third 
root is easy to find 



^ ^(15p2 - 2Sopo T 4:k^/9n^ - G^oPoP^) 
ax = = — , = . (17) 

2v^V9fcp(p2 + 2^oPo) ± v^(3/i2 _ 2£^p,)l 
It corresponds to the radius where the universe turns around or bounces. 
Now we divide our discussion into two cases: e = +1, and e = — 1. 
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A. e = +l 



Since the oscillating universe exists only in the case of e = +1 and k = 1, we first focus 
on this case. 



1. k = l 

By introducing d = \/ gm(^ — f^'^^o, Eq- ( 1TT|) becomes 

Apparently, when \^ = 0, we get a cubic equation of a, which can be expressed as 

~ "^^"I^T^TT^^FT*^^ ~ «mm)(a - amax){a -dx) = , (19) 
ofi yd + /i j 

with dmin, CLmax and fir being three solutions. Assuming < a^m < S^ax < 5^, if 
V(a) < in a G [amm,amax] and the equality holds when d = dmin and d = dmax, the 
universe oscillates between two finite scales; if V{d) < in a G [dx, oo) and the equality 
holds when a = ot, it corresponds to a bounce scenario and the universe bounces at dx- 
For a simple example, let g'^ = Sfi^So (the boundary T (Eq. f|T3l) ) in Fig. ([2]) for obtaining 
the oscillating universe) in the potential (Eq. (fTSll ) , we have 

(^min = CtO ) (20) 



(21) 



ot = 7= ) • (22) 

For the case < |, we have < amin < CLmax < o-t- In Fig. ([T]), we plot the evolutionary 
curve of V{a). It is easy to see that V{a) < in a G [a^jn, amax] and a G [a^-, oo), which 
means there is an oscillating universe (a G [a^jn, Omai]) or a bouncing one (a G [ax, oo)). 
Thus, if the universe is in the region [amin, cimax] initially, it may undergo an oscillation. 
After a number of oscillations, it may evolve to the bounce point through quantum 



FIG. 1: The potential V{a) for a matter dominated universe under conditions = Sfj.'^So and 
< |. The constants are set as ^ = 1, po = 0.1, gm = 5.61249, and £q = 10.5. The radii 
of the oscihation are Omm = 1.87083 and amax = 2.97382. The period of this oscillation is 
T = 18.3324. The largest root ut = 5.04402 corresponds to a turning radius of a bouncing 
universe. 

tunneling. If the universe contracts initially from an infinite scale, it can turn around at 
qt and then expand forever. 

For a general case, we find that there is an oscillating universe if the following conditions 
are satisfied 

9 

o 

9 3 

-P^ < £oPo < , g;n<9m< 9^ (23) 

where is defined in Eq. (fT5|) . Using above equations, we obtain the allowed region in 
(^imp.^ ^T^) plane (Fig. (|2])) for an oscillating universe or a bouncing one. The boundaries 
curves T± are defined as gm = gti (Eq. ( fTSil ) and curve V a.s gm = y/SSofi'^ (Eq. (fT3ll). 
The period of an oscillation can be calculated through 

Z""™- da , 



mm 



where amax and are the maximum and minimum radius of the oscillating universe. 

In Fig. ([3]), we give the evolutionary curve of the potential with the model parameters 
satisfying Eq. (!23|) . From this figure, we find that there is an oscillating universe between 
amin and amax, or a bouncing one in [ar,oo). Therefore, a similar cosmic evolution as 
shown in Fig. ([1]) is obtained. 
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FIG. 2: Phase diagram of spacetimes in ( , -^) plane for a matter dominated universe. An 
oscillating universe is found in Region 1. A bounce one is found in Regions 1, 2, and 3. The 
unstable and stable static universes exist on the boundaries and r_, respectively. 




FIG. 3: The potential V{a) for an oscillating universe or a bouncing one in Region 1 of Fig. ([2]). 
The constants are chosen as = 1, po = 0.1, gm = 5.75, and £q = 11.25. The radii are 
Omin = 1.96042, Umax = 3.11647 and ax = 4.70572. The period of an oscillation is T = 21.9044. 

On the r+ curve, the unstable static universe appears. The solution is 

as = , (25) 

with a+ given in Eq. ffT6|) . which is a double solution of V{a) = 0. If < Sopo < 
the third root is 

am = , (26) 

with given in Eq. f lT7|) . In Fig. @ we plot the evolution of the potential. When a = as, 
both V and V vanish, but this as solution is unstable. Therefore, in the < £opo < 
case , the universe can oscillate between and as, and it can also evolve directly from a^, 
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to oo or evolve to oo after some oscillations with no need of quantum tunneling to make 
it happen. If the universe contracts initially from an infinite scale, it can turn around 
at as, or pass through as and bounce at a^, then oscillate between and as- For the 
£oPo < fyU^ case, if the universe initially evolves from oq, it can expand to as, and then 
it can further expand to infinity or turn around. Once the universe bounces at as and 
contracts to oq, there will appear a new singularity since at Oq the Hubble parameter 
H{ao) is nonzero. That is, when the universe contracts to ag and then expands, it has to 
evolve discontinuously at oq. 




FIG. 4: The potential V{a) in a matter dominated universe for an unstable static universe (on 
boundary r_|_ of Fig. ([2])). The constants are chosen /x = 1, po = 0.1, So = 11-5, Qm = 5.84184 
(left), and fi = 1, po = 0.1, Sq = 11, gm = 5.77729 (right). The radii are am = 1.96954, 
as = 3.85103 (left) and oq = 1.90400, as = 3.89410 (right). 

A stable static universe or a bouncing one exists on the r_ curve. The stable static 
solution is 

as = ag (27) 

with aj given in Eq. ( 1T6|) . while the turning radius ar of a bouncing universe is given by 

ar = a+ (28) 

with a J given in Eq. f|T7|) . Fig. ([5]) gives the potential with model parameter in the F_ 
curve. There are two solutions {as, ot) for ^ = 0. At a = as, both V and V are equal 
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to zero and apparently a = as corresponds to a stable solution. Thus, the universe can 
stay at this finite radius past-eternally. When a > ax, V < 0, which corresponds to a 
bouncing universe. If the universe stays at initially, after a long time, it can quantum 
mechanically tunnel to the bounce point ot and then expand. If the universe contracts 
initially from an infinite scale, it will turn around at ar- 

tv 

0.05 



-0.05 



-0.10 




FIG. 5: The potential V{a) in a matter dominated universe for a stable static universe and 
a bouncing one (the (r_) curve of Fig. ([2]))- The constants are chosen as fj, = 1, po = 0.1, 
dm = 5.75374, and Sq = 12. The radius of stable static universe is as = 2.35114 and the bounce 
radius is ut = 5.20431. 



Fig. ([6]) shows the evolution of the potential V{a) with the model parameters in Region 

2 of Fig. ([2]). From this figure, we find that a bouncing universe is obtained since ^ < 
in a G [aT,oo). In addition, V <0 in a G [aQ,am], but = occurs only at a = a^- 
Thus, if the universe turns around at and contracts to ag, as shown in the left panel 
of Fig. (jl]), there appears a singularity at oq. Of course, if the universe evolves from ao to 
ttm and then quantum tunnels to directly, this singularity can be avoided. 

Fig- © shows the evolution of the potential V{a) with the model parameters in Region 

3 of Fig. (j2]). Apparently, a bouncing universe is obtained. 

We plot Fig. ([HD to give the evolution of the potential with the model parameters in 
Region 4 of Fig. ([2]). In this case < for a > Cq. Thus an FS ^ oo type is obtained. 



11 



FIG. 6: The potential V{a) in a matter dominated universe with model parameters in Region 2 
of Fig. ([2]) . The constants are chosen as /i = 1, />o = 0.1, gm = 4.47214, and So = 5. The radii 
are ao = 1.11803, = 2.23607, and ar = 7.63441. 
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FIG. 7: The potential V{a) in a matter dominated universe for a bouncing universe (Region 
3 of Fig. ([2])). The constants are /i = 1, po = 0.1, Qm = 3.16228, and Sq = 10. The radius is 
OT = 15.5259. 
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FIG. 8: The potential V{a) in a matter dominated universe for model parameters in Region 4 
of Fig. ([2|). The constants are chosen as n = 1, po = 0.1, Qm = 7.07107, and So = 10. 
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2. k = -l 



In this case, the potential is given by 

1 9n 



V{a) 



^0 



(29) 



3/2^ a 3/2^ 

Fig. shows the evolution of this potential, from which one can see that, as Fig. ([H]), 
the potential V{a) is always negative for a > Oq, which means that the cosmic evolution 
type is FS ^ oo. 




10 12 14 



FIG. 9: The potential V{a) in a matter dominated universe for the case e = +1, and k = — 1. 
The spacetime will always expand. The constants are chosen as n = 1, po = 0.1, gm = 10 and 
£o = 10. 



B. e = -l 

In this case, we get 

Since V'{a) is always positive, the potential is an increasing function. Therefore, the 
cosmic evolution will be simpler than the case of e = +1. 



1. k = l 

From Fig. (fTOl) one can see there is no solution for V < if Qm < y/SSofJ^'^- When 
gm > y/SSofj,"^, if tto < a < am, < 0, while = occurs only at am- So, the universe 

13 



(a) g^>^ 3 So /^^ (b) 8m<yl 3 So /^^ 




FIG. 10: The potential l^(a) in a matter dominated universe for e = — 1 and = 1 with 

gm > \/3<5o/i^ (left) and < V^'^oA*^ (right). The constants are chosen as /i = 1, = 0.1, 
So = 10, = 12 (left) and ^ = 1, = 0.1, Sq = 10, Qm = 5 (right). 

can evolve between oq and am, but there is a singularity at oq since H{ao) 7^ at this 
point, which means that this cosmic evolution type is FS =^ FS. When Qm = v^3^/i^, 
there is only one point a = Oq for < 0. 

2. k = -l 
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FIG. 11: The potential V(a) in a matter dominated universe for the case of e = —1 and k = — 1. 
The constants are e = —1, k = —1, 7/ = 1, /i = 1, /jq = 0.1, £0 = 10, and Qm = 15. 

From Fig. ( fTTl) we find that, when cto < < c^m, < 0. So, a FS* ^ FS* type is 
obtained. 



14 



IV. THE EVOLUTION OF A RADIATION-DOMINATED UNIVERSE IN THE 
DGP BRANEWORLD 

In this section, we discuss the case where the universe is dominated by radiation 
{to = I). Thus, the cosmic energy density can be expressed as p = ^i-, and the potential 
becomes 

na)={k-^,M-,^e,,))-^,'^. (31) 

It is easy to see that Qr > p^i^o is required. As in the previous section, we divide our 
discussions into two cases: e = +1, and e = — 1. 

A. e = +l 

1. k = l 




FIG. 12: The potential V{a) in a radiation dominated universe for e = 1 and k = 1 with 
9r > 1^ + fJ-'^So (left) and gr < ^ + Ai^<?o (right). The constants are chosen as /x = 1, = 0.5, 
£q = 5, Qr = 20 (left) and fj, = 1, po = 0.5, So = 5, gr = 8 (right). The radius where the universe 
turns around is ax = 2.51185. 

As shown in Fig. fll2p . when > |^ + /i^i^o; the potential is always negative and the 
type of the cosmic evolution is BB ^ oo or oo ^ BC. While, for /i^^o < 9r < f^ + p^^o, 
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the potential will turn to be positive from negative at the radius 



aT=J , (32) 

3/i2 - ^/2po{gr - fi'^So) 



Thus, the cosmic evolution type is BB =^ BC . 



2. k = -l 



V 



2 4 s a 




FIG. 13: The potential V{a) in a radiation dominated universe for e = 1 and k = —1. The 
constants are chosen as ^ = 1, po = 0.5, Qr = 5, and Sq = 5. 

From Fig. f lT3|) we can see that in this case the potential is always negative and the 
cosmic evolution type is BB ^ oo or oo ^ BC. 



B. e = -l 



1. k = l 



The type of the cosmic evolution now is BB =^ BC as can be seen from Fig. (IHj). The 
turning point is 



aT = J , (33) 

3/x2 + ^/2po{gr - fi^So) 
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FIG. 14: The potential V{a) in a radiation dominated universe for e = —1 and k = 1. The 
constants are chosen as /i = 1, po = 0.5, gr = 5, and So = 5. The bounce radius is ar = 1.29099. 



(a) gr>fi^So+ 



2po 



(b) gr<I^^So+ 




2po 



3. 



FIG. 15: The potential V{a) in a radiation dominated universe for e = — 1 and k = —\ with 



(7r. > 2^ + £-Q (left) and < 2^ + "^^o (right). The constants are chosen as // = 1, = 0.5, 
^0 = 5, 5r = 20 (left) and ^ = 1, po = 0.5, <fo = 5, = 8 (right). The bounce radius is 
ar = 3.87298. 



2. k 



-1 



From Fig. ( !T5|) we obtain that when (yf^. > -j^^ — h /i^£^o the cosmic evolution type is 



2po 

55 =^ BC. The radius where the universe turns to contract is 



-3/^2 + ^j2po{gr - fi^So 



(34) 



When fi^So < gr < ^ + /i^^^o, the cosmic evolution type is BB ^ 00 or 00 ^ BC since 



the potential is always negative. 



17 



V. THE EVOLUTION OF A VACUUM-DOMINATED UNIVERSE IN THE 
DGP BRANEWORLD 



If the universe is dominated by vacuum energy [u = —1), the cosmic energy density is 
a constant. We denote it hj p = g^, and the potential becomes 

V{a) = k - jj^dvO^ - -^V^PoiOva^ - P^^o) ■ (35) 



9v 

e = — 1 and 



Clearly, a > = .yi^^ is needed, and usually V{ao) ^ except for the case k = 1, 



9v = ■ (36) 

•^0 



A. e = +l 



1. k = l 




FIG. 16: The potential V{a) in a vacuum dominated universe for e = 1 and k = 1 with > 
(left) and < ^ (right). The constants are chosen as /i = 1, po = 0.1, So = 10, = 15 (left) 
and fi = 1, po = 0.1, So = 3, gv = 2 (right). And the radius where the universe turns around is 
ar = 1.1547. 

As plotted in Fig. f|T6|) . the potential V{a) is a decreasing function of a. If > the 
potential V{a > ao) is always negative and the type of the cosmic evolution is FS =^ oo. 
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If gv < Viar) = 0, thus, we get a bouncine; universe. The bounce radius is 



Or 



(37) 



2. k = -l 



We find from Fig. ( IT71) that the cosmic evolution type is FS =^ oo. 




FIG. 17: The potential V{a) in a vacuum dominated universe for e = 1 and A; = —1. The 
constants are chosen as /i = 1, = 0.1, gv = 15, and £q = 10. 



B. e = -l 



i. /t = 1 



In this case, there exists a static universe if g^ and other parameters satisfy a relation 



9v 



So 



2po = 



(3^ 



which is obtained by combining V{a) = and V'{a) = 0. Using the above equation, one 
can obtain a static state solution 



as 



(39) 



Using Eqs. fl36| [391) . we can depict all cosmic evolution types in the (^, ^) plane, 
which is shown in Fig. (fT8|) . On the green line in Fig. (fT8|) . which is determined by Eq. (136|1 . 
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FIG. 18: Phase diagram of spacetimes in -^) plane for a vacuum dominated universe with 
e = — 1 and k = 1. A. bouncing universe is found in Regions 1, and 2. oq =^ oo type of universe 
are found in Region 3. The unstable and stable static universes exist on the curves Tsu and F, 
respectively. 

we can get a stable static universe with as = ao and a bouncing one as shown in Fig. flT9|) . 
Thus, the universe can originate from a stable Einstein static state, which means that 
the universe stays at this stable state past-eternally and then enters a expanding phase 
through quantum tunneling. If the universe contracts initially from an infinite scale, it 
will bounce at qt- The radius of the stable static universe is 



While, on the blue line of Fig. ( !T8|) . which is determined by Eq. ( 138|) . we obtain an 
unstable static universe. Its radius is 



We plot the effective potential in Fig. fl20|) . From which, one can see that the cosmic 
evolution type is similar to that shown in the right panel of Fig. 

Fig. ( l2Ti) shows the evolution of the potential V{a) with the model parameters in 
Region 1 of Fig. (ITSl) . We find that a bouncing universe is obtained and the bounce 




(40) 



and the bounce radius is 




(41) 




(42) 
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FIG. 19: The potential V{a) for a stable static universe or a bouncing one in a vacuum dominated 
phase (the green line of Fig. ([18]) ). The constants are chosen as fi = 1, po = 0.05, = 0.9, and 
So = 10. The radii are as = 1.82574 and ot = 2.04124. 
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FIG. 20: The potential V{a) for an unstable static universe in a vacuum dominated phase (the 
blue line of Fig. p^ ). The constants are chosen as fi = 1, po = 0.1, (7^ = 1.1, and £0 = 10. The 
radii are oq = 1.73640 and as = 1.82574. 




FIG. 21: The potential V{a) in a vacuum dominated universe with the model parameters in 
Region 1 of Fig. (fTHj) . The constants are chosen as ^ = 1, po = 0.1, gv = 0.5. The bounce radius 
is ar = 3.8941. 
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radius is 



gvi9v - '^Po) 



(43) 
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FIG. 22: The potential V{a) in a vacuum dominated universe with the model parameters in 
Region 2 of Fig. ([TH]) . The constants are chosen as /i = 1, po = 0.3, = 1.4, and £q = 10. The 
radii are oq = 1.63481, am = 1.71222, and ax = 2.13736. 

In Fig. fl22|) . we plot the potential V{a) with the model parameters in Region 2 of 
Fig. (fTSl) . A similar result as shown in Fig. (jH]) is obtained. The expressions for and 

ax are 



' 35(^/^2 _ y/2g^pQp^{9p^ - Sogy + 2SoPo) 



g^i^g^ - 2po) 



(44) 



ax 



1 3gyp^ + ^/2gypop'^{9p^ - Sogv + 2£oPo) 



(45) 



9vi9v - '2po) 

Fig. fl23!) shows the evolution of the potential V{a) with the model parameters in 
Region 3 of Fig. (fT8|) . which corresponds to the cosmic evolution type: FS =^ oo. 



2. k = -l 



It is easy to see from Fig. fl24|) that V{a) is always negative. So the cosmic evolution 
type is FS =^ oo. 
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FIG. 23: The potential V{a) in a vacuum dominated universe with the model parameters in 
Region 3 of Fig. ([18]). The constants are chosen as = 1, po = 0.1, gv = 15, and £o = 10. 
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FIG. 24: The potential V{a) in a vacuum dominated universe with e = —1 and k = —1. The 
constants are chosen as fj, = 1, po = 0.1, = 10, and Sq = 10. 

VI. CONCLUSIONS 

In this paper, we have studied all possible cosmic evolutions in the DGP braneworld 
scenario with a method in which the dynamics of the scale factor is treated like that 
of a particle in a "potential". The effect of the dark radiation on the cosmic evolution 
is considered. By assuming that the cosmic energy component is pressureless matter, 
radiation or vacuum energy, respectively, we find that, in the matter or vacuum energy 
dominated case, the universe does not originate from the big bang singularity and its scale 
factor has a minimum value oq. Thus the classical singularity problem can be avoided. 
However, there may appear a new singularity at in the sense that when the universe 
bounces or contracts to this point and then expands, its evolution will be discontinuous 
as H{ao) 7^ 0. However, in some circumstances, there exists a stable Einstein static state 
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or a bouncing universe to avoid the new and classical singularities. If the universe is in 
the Einstein static state initially, it can stay there past-eternally and evolve to the bounce 
point through quantum tunneling. If the universe contracts initially from an infinite scale, 
it can turn around at the bounce point and then expand forever. Therefore the cosmic 
evolution is nonsingular. In addition, in the matter dominated case, there also exists an 
oscillating universe to avoid the singularity problem as long as the model parameters are 
in some specific regions (shown in Fig. (|2])). If the cosmic scale factor is in the oscillation 
region initially, the universe may undergo an oscillation. After a number of oscillations, it 
may evolve to the bounce point through quantum tunneling. In the radiation dominated 
case, the universe may originate from the big bang singularity, but a bouncing universe 
which avoids this singularity is also possible. 
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